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PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Teehnieal  Note  (CHETN)  provides  empirieal 
equations  to  estimate  flow  velocity  and  flow  thickness  resulting  from  irregular  waves  overtopping  a 
trapezoidal-shaped  earthen  levee.  The  CHETN  summarizes  published  European  research  with  the 
goal  of  providing  practicing  engineers  with  design  guidanee  in  a  eoneise  and  understandable  form. 
Worked  examples  illustrate  application  of  the  empirical  equations. 

INTRODUCTION  AND  BACKGROUND:  Erosion  protection  (armoring)  systems  placed  on  the 
crown  and  protected-side  slope  must  resist  the  forces  of  fast-flowing,  turbulent  water  that  has 
overtopped  the  levee  erest.  Eigure  1  illustrates  the  three  overtopping  cases  that  might  occur: 
(a)  overtopping  by  wind-generated  waves  when  the  still  water  level  is  beneath  the  elevation  of  the 
crest,  (b)  overflow  by  water  levels  above  the  levee  crest,  but  without  wave  activity,  and 
(c)  overtopping  by  combined  waves  and  storm  surge. 


a)  Wave  Overtopping 


b)  Surge  Overtopping 


Figure  1 .  Overtopping  scenarios  for  earthen  levees. 
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The  guidance  presented  in  this  CHETN  pertains  only  to  irregular  wave  overtopping  while  the  surge 
level  (still  water  elevation)  is  lower  than  the  levee  crest  elevation,  as  shown  in  the  top  sketch  of 
Figure  1 .  The  main  difference  between  wave  overtopping  and  steady  flow  overtopping  (shown  in  the 
middle  sketch  of  Figure  1)  is  the  unsteady  flow  that  characterizes  periodic  wave  overtopping.  As 
each  wave  overtops,  it  has  a  velocity  across  the  levee  crest  that  can  exceed  the  velocity  of 
comparable  surge  overtopping.  Thus,  unprotected  soil  on  the  levee  crest  that  is  stable  for  low  levels 
of  surge  overtopping  may  erode  when  waves  overtop. 

For  typical  slopes  on  the  protected  side  of  levees,  overtopping  waves  can  create  critical  flow 
conditions  near  the  leeward  edge  of  the  levee  crest,  resulting  in  supercritical  wave  overtopping  flow 
on  the  backside  slope.  However,  this  flow  condition  is  unsteady,  and  peak  velocities  are  sustained 
for  only  a  brief  time.  In  addition,  the  unsteady  discharge  over  the  crest  results  in  a  limited 
overtopping  volume.  Consequently,  any  erosion  on  the  backside  slope  due  to  wave  overtopping  is 
intermittent,  and  the  erosion  rate  will  vary  with  overtopping  intensity. 

SUMMARY  OF  EUROPEAN  EXPERIMENTS:  Experiments  have  been  conducted  in  Europe  at 
small  and  large  scale  with  the  aim  of  quantifying  the  wave  overtopping  flow  parameters  on  the  inner 
slope  of  dikes  and  levees  (Schiittrumpf  et  al.  2002;  van  Gent  2002;  Schiittrumpf  and  van  Gent  2003; 
Schiittrumpf  and  Oumeraci  2005).  The  European  researchers  developed  analytical  expressions  to 
represent  the  velocity  and  flow  depths  at  the  edge  of  the  crest  on  the  flood  side  {ua,  Ha),  at  the  edge 
of  the  crest  on  the  protected  side  (ub,  hs),  and  down  the  backside  slope  {ush,  hsb),  as  illustrated  in 
Figure  2. 


Figure  2.  Wave  overtopping  definition  sketch  (after  Schuttrumpf  and  Oumeraci  2005). 


The  key  parameters  necessary  for  estimating  the  flow  velocities  and  depths  are  the  levee  freeboard, 
Rc,  the  runup  elevation  exceeded  by  2  percent  of  the  waves,  Ru2%,  and  a  friction  factor,  fp,  that 
accounts  for  frictional  energy  loss  as  the  overtopping  wave  travels  across  the  levee  crest  and  down 
the  protected-side  slope. 

Independent  laboratory  experiments  were  conducted  in  The  Netherlands  (van  Gent  2002)  and  in 
Germany  (Schuttrumpf  et  al.  2002).  These  two  studies  produced  similar  estimation  analysis 
techniques,  with  only  minor  differences  in  the  details.  A  joint  paper  (Schuttrumpf  and  van  Gent 
2003)  reconciled  the  differences  to  the  extent  possible. 
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Van  Gent’s  (2002)  small-scale  experiments  had  a  1 ;  100  foreshore  slope  with  a  1 :4  slope  on  the  flood 
side  of  the  dike.  Two  levee  crest  widths  (0.2  and  1.1m)  were  tested  with  two  protected-side  slopes 
(1:2.5  and  1:4)  to  give  four  dike  geometries  using  a  smooth  dike  surface.  A  fifth  test  series  was 
conducted  with  a  rough  surface.  Velocity  and  flow  thickness  were  measured  at  the  flood-side  and 
protected-side  edges  of  the  crest  and  at  three  locations  spaced  down  the  protected-side  slope.  Micro¬ 
impellers  were  used  to  measure  velocity.  Eighteen  irregular  wave  tests  were  performed  for  the 
different  dike  geometries,  ten  with  single-peaked  spectra  and  eight  with  double-peaked  spectra. 
Incident  wave  conditions  were  determined  by  measuring  the  generated  waves  without  the  structure 
in  place,  and  applying  the  Mansard  and  Funke  (1980)  frequency-domain  method  to  remove 
reflection  caused  by  the  dissipating  beach  profile.  Van  Gent  (2002)  used  the  wave  parameter  //1/3  in 
the  analysis,  but  did  not  indicate  how  this  time-domain  parameter  was  determined  from  the 
frequency-domain  value  of  Hmo  found  from  the  reflection  analysis.  Wave  period  was  specified  as 
mean  period  I'm- 1.0,  and  it  was  estimated  from  the  moments  of  the  incident  wave  frequency  spectra. 
The  mean  period  is  reported  to  better  represent  double-peaked  spectra. 

The  experiments  by  Schiittrumpf  et  al.  (2002)  included  both  small-  and  large-scale  tests.  The  small- 
scale  tests  used  three  flood-side  slopes  (1:3, 1:4,  and  1:6),  a  crest  width  of  0.3  m,  and  five  protected- 
side  slopes  (1:2,  1:3, 1:4, 1:5,  and  1:6).  A  total  of 270  small-scale  tests  were  run  using  regular  waves 
and  irregular  waves  conforming  to  the  JONSWAP  spectrum.  Flow  depths  were  measured  with 
resistance  wave  gauges,  and  overtopping  flow  velocity  was  recorded  using  micro-impellers.  For  the 
large-scale  tests  the  flood-side  slope  was  1 :6,  the  crest  width  was  2  m,  and  the  protected-side  slope 
was  1:3.  A  total  of  250  large-scale  model  tests  were  run  using  some  regular  waves,  but  mostly 
irregular  waves.  Flow  depth  and  velocity  were  measured  using  wave  gauges  and  micro-impellers. 

Wave  data  from  tests  by  Schiittrumpf  et  al.  (2002)  were  analyzed  in  the  frequency  domain  using  the 
reflection  method  of  Mansard  and  Funke  ( 1 980).  The  time-domain  wave  height  parameter  Hy^  was 
used  in  their  overtopping  analysis  with  the  conversion  from  the  frequency  domain  wave  height  given 
as  Hi/3  =  0.94  H^q  (Schiittrumpf  2006,  personal  communication).  This  conversion  is  a  little  odd 
because  we  should  expect  Hy^  to  be  greater  than  H^o  for  shallow  water  waves.  Also,  the  conversion 
is  valid  strictly  for  these  tests  and  not  in  general  because  it  was  determined  for  wave  flume  data  with 
a  constant  water  depth  for  all  tests.  The  wave  period  was  specified  as  the  mean  wave  period,  and  it 
was  determined  from  the  calculated  incident  wave  spectra  by  the  simple  relationship  !„  =  0.88  Tp 
(Schiittrumpf  2006,  personal  communication). 

FLOW  PARAMETERS  AT  THE  FLOOD-SIDE  LEVEE  CREST  EDGE:  Flow  depth  of  the 
runup  on  the  flood-side  slope  was  assumed  to  be  a  linear  decrease  from  the  still  water  level  to  the 
highest  elevation  of  runup.  At  the  flood-side  edge  of  the  levee  crest  (denoted  by  the  subscript  letter 
A  in  Figure  2),  the  flow  parameters  are  given  by  the  equations 


"■A2%  _ 
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{^u2%  ^c) 
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and 


^A2% 

4^. 


=  C 


Au2%^ 


{^u2%  ^c) 

H, 


where 


hA2% 

Hs 

CAh2% 

Ru2% 

Rc 

Ua2% 

g 

Cau2% 


peak  flow  depth  exeeeded  by  2%  of  the  waves 

signifieant  wave  height  [~  Hn,o\ 

empirieal  depth  eoeffieient  determined  from  test  data 

runup  elevation  exceeded  by  2%  of  the  waves 

crest  freeboard  [=  crest  elevation  minus  still  water  elevation] 

flow  depth-averaged  peak  velocity  exceeded  by  2%  of  the  waves 

acceleration  of  gravity 

empirical  velocity  coefficient  determined  from  test  data 


(2) 


The  values  of  hA2%  and  ua2%  were  determined  from  the  peaks  of  the  overtopping  wave  time  series, 
and  these  parameters  represent  the  levels  exceeded  by  only  2  percent  of  the  total  waves  during  the 
tests.  For  example,  if  a  test  had  1000  waves,  perhaps  only  200  waves  overtopped  the  crest.  The 
2  percent  exceedance  level  would  be  the  level  exceeded  by  20  of  the  1000  waves  (0.02  x  1000),  but 
this  is  10  percent  of  the  overtopping  waves.  Schiittrumpf  et  al.  (2002)  also  provided  coefficients  for 
the  average  overtopping  parameters  /z^50%  and  ua50%-  All  of  the  equations  pertain  to  the  maximum 
velocity  at  the  leading  front  of  the  overtopping  wave.  Flow  velocities  and  depths  associated  with  a 
single  wave  decrease  after  passage  of  the  wave  front. 

Note  in  Equations  1  and  2  that  significant  wave  height  Hs  in  the  denominator  cancels  on  both  sides 
of  the  equations.  Thus,  the  flow  depth  is  directly  proportional  to  the  difference  between  the  2-percent 
runup  and  levee  freeboard,  and  the  depth-averaged  flow  velocity  is  proportional  to  the  square  root  of 
this  difference.  Wave  parameters  enter  into  the  estimation  of  flow  depth  and  velocity  at  the  flood- 
side  crest  edge  through  the  estimation  of  the  2-percent  runup  parameter  As  noted  by  van  Gent 
(2002),  the  calculated  Ru2%  is  a  fictitious  value  in  cases  where  runup  exceeds  the  structure  freeboard. 
It  is  the  level  that  would  be  exceeded  by  2  percent  of  the  waves  if  the  front  slope  were  continued 
upward  indefinitely. 

The  values  of  the  empirical  coefficients  determined  for  the  two  studies  are  given  in  Table  1.  The 
superscripts  following  each  number  refer  to  the  references  given  in  the  footnotes  to  Table  1. 
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Table  1.  Empirical  coefficients  for  flood-side  crest  edge  flow 
parameters. 

Coefficient 

Schuttrumpf 

van  Gent 

C/t/72% 

0.33  and  0.22 

0.15 

1.55^  and  1.37^ 

1.30 

C/t/750% 

0.17^'^ 

— 

C4u50% 

0.94 

— 

'  van  Gent  (2002). 

^  Schuttrumpf  et  al.  (2002). 

^  Schuttrumpf  and  van  Gent  (2003). 

Schuttrumpf  and  Oumeraci  (2005). 

The  coefficient  C^/,2%  is  a  constant  that  is  actually  equal  to  a  slope-dependent  constant,  C2,  divided 
by  tan  0,  where  0  is  the  flood-side  structure  slope.  Values  of  C2  given  in  the  various  papers  are  used 
in  an  equation  slightly  different  from  Equation  1 .  The  value  for  C^/,2%  given  by  Schiittrumpf  was 
revised  from  0.33  to  0.22  in  the  most  recent  paper  (Schiittrumpf  and  Oumeraci  2005),  and  this 
probably  represents  a  better  value — as  shown  by  the  data  plot  given  in  their  paper  and  the  fact  that  it 
is  closer  to  the  value  obtained  by  van  Gent.  The  value  of  Cau2%  =  1.55  is  derived  from  a  table  in 
Schiittrumpf  et  al.  (2002)  that  associated  this  coefficient  with  large-scale  tests.  A  coefficient 
associated  with  the  10-percent  exceedance  level  can  also  be  derived  from  the  same  table  as  Cau\o% 
=  1.37  for  large-scale  tests.  In  Schiittrumpf  and  van  Gent  (2003),  the  value  of  CauIX  =  1-37  was 
reported,  and  this  is  thought  to  be  a  typographical  error.  The  correct  value  should  have  been  Cau2% 
=  1.55. 

Schiittrumpf  and  van  Gent  (2003)  attributed  differences  in  empirical  coefficients  to  different  dike 
geometries  and  instruments,  but  noted  the  differences  are  not  too  great.  Van  der  Meer  et  al.  (2006) 
suggested  an  error  in  measurement  or  analysis  might  have  caused  the  factor-of-two  difference  seen 
for  the  coefficient  C^/,2%,  but  the  revised  value  of  C^/,2%  =  0.22  brings  the  results  closer.  A  more 
probable  cause  for  variation  might  be  in  the  method  each  investigator  used  to  estimate  the  value  of 
2-percent  runup,  Ru2%-  Van  Gent  (2002)  estimated  Ru2%  using  a  formula  he  developed  earlier  (van 
Gent  200 1 ),  which  uses  //1/3  and  Tm-o.  1  as  the  wave  parameters.  Schiittrumpf  estimated  Ru2%  using  the 
equations  of  de  Waal  and  van  der  Meer  ( 1 992),  with  wave  height  and  wave  period  Tm  instead  of 
spectral  peak  period  Tp.  Both  formulas  give  reasonable  estimates  that  fall  within  the  scatter  of  the 
2-percent  runup  data;  thus,  whichever  formula  is  selected  for  calculating  Ru2%,  the  estimates  for 
overtopping  flow  parameters  should  be  similar. 

Until  further  clarification  becomes  available,  it  is  recommended  that  values  of  CAhiy«  =  0.22  and 
Cau2»/o  =  1.55  be  used  in  Equations  1  and  2  to  estimate  the  overtopping  flow  parameters  associated 
with  the  flow  depth  and  velocity  exceeded  by  2  percent  of  the  incoming  waves. 

FLOW  PARAMETERS  AT  THE  PROTECTED-SIDE  LEVEE  CREST  EDGE:  Overtopping 
waves  flowing  across  the  dike  or  levee  crest  decrease  in  height,  and  the  velocity  decreases  as  a 
function  of  the  surface  friction  factor,^)’.  Elow  depth  (or  thickness)  can  be  estimated  at  any  location 
on  the  crest  with  the  equation 
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^  X  ^ 

-Q- 

'  B 


(3) 


where  C3  is  an  empirieal  eoeffieient,  Xc  is  distance  along  the  crest  from  the  flood-side  edge,  and  5  is 
the  horizontal  crest  width.  The  flow  thickness  at  the  protected-side  crest  edge  (denoted  by  the 
subscript  letter  B  in  Figure  2)  is  given  whenxc  =  B.  Different  values  of  the  coefficient  were  given  in 
the  various  publications,  i.e.,  based  on  the  2-percent  exceedance  levels  C3  =  0.89  for  TMA  spectra 
and  C3  =  1 . 1 1  for  natural  spectra  (Schiittrumpf  et  al.  2002);  C3  =  0.40  and  0.89  (Schuttrumpf  and  van 
Gent  2003);  and  C3  =  0.75  for  irregular  and  regular  waves  (Schuttrumpf  and  Oumeraci  2005).  The 
factor-of-two  difference  between  van  Gent  and  Schuttrumpf  was  attributed  to  the  difference  in 
estimating  wave  runup. 

For  levee  calculations,  it  is  recommended  that  a  value  of  C3  =  0,75  be  used  in  Equation  3  on  the 
assumption  that  earlier  values  have  been  superseded  by  publication  of  the  2005  journal  article.  The 
2-percent  runup  elevation  should  be  estimated  using  the  runup  formulas  of  de  Waal  and  van  der 
Meer  (1992)  or  Hughes  (2004).  If  van  Gent’s  (2001)  method  for  estimating  wave  runup  is  used,  it 
would  be  more  appropriate  to  use  a  value  of  C3  =  0.40.  Note  that  Equation  3  is  applicable  for 
estimating  /zs50%  if  the  flow  depth  hA50%  is  used  instead  of  hA2%-  In  fact,  Schuttrumpf  and  Oumeraci 
(2005)  presented  only  the  50-percent  exceedance  values. 

Elow  velocity  along  the  dike  crest  exceeded  by  2  percent  of  the  waves  is  given  by  a  similar  equation; 


^B2% 


f 

M^20/.exp 

V 


fp 

2  hg2%  j 


(4) 


where  fp  is  the  banning  friction  factor  appropriate  for  the  levee  crest  surface,  and  hsix  is  the  flow 
depth  at  that  location  on  the  crest  obtained  via  Equation  3.  At  the  protected-side  crest  edge,  evaluate 
Equation  4  with  Xc  =  5.  Van  Gent  (2002)  had  a  different  expression  for  ub2%',  however,  in 
Schuttrumpf  and  van  Gent  (2003),  both  authors  agreed  on  Equation  4.  A  theoretical  derivation  for 
Equation  4  was  given  in  Schuttrumpf  and  Oumeraci  (2005). 

ESTIMATION  OF  FRICTION  FACTOR:  The  Fanning  friction  factor  has  a  significant  influence 
on  flow  velocity  across  the  crest  and  down  the  protected-side  slope.  The  small-scale  experiments  of 
Schuttrumpf  et  al.  (2002)  had  a  structure  surface  constructed  of  wood  fiberboard,  and  the  friction 
factor  on  the  crest  was  determined  experimentally  to  be^F  =  0.0058  (Schuttrumpf  and  Oumeraci 
2005).  The  structure  in  the  companion  large-scale  experiments  was  constructed  with  a  bare, 
compacted  clay  surface,  and  experimental  results  gave  the  friction  factor  as  fF  =  0.01  (Schuttrumpf 
et  al.  2002).  Schuttrumpf  and  Oumeraci  (2005)  also  list  the  representative  value  for  friction  factor  on 
the  protected-side  slope  as^F  =  0 . 02  (smooth  slopes),  and  Cornett  and  Mansard  ( 1 994)  give ^=0.1 
-  0.6  (rough  revetments  and  rubble-mound  slopes).  Grass-covered  slopes  would  probably  have  a 
friction  coefficient  above  fF=  0.01. 

Determination  of  an  appropriate  value  of  friction  factor  for  various  grass  coverings  and  armoring 
alternatives  may  be  difficult  because  of  the  lack  of  published  values.  As  a  first  approximation,  an 
estimate  can  be  made  if  a  representative  value  of  Manning’s  n  is  known  for  a  particular  slope  surface 


6 


ERDC/CHLCHETN-lll-77 
January  2008 


or  armoring  product.  Manning’s  n  can  be  related  to  the  Chezy  eoeffieient,  Cz,  by  the  expression  (e.g., 
Henderson  1966) 


C,  = 


R 


1/6 


(5) 


where  R  is  the  hydraulie  radius,  and  Manning’s  n  is  given  in  metrie  units.  For  wide  ehannels,  R  is 
essentially  the  same  as  the  depth,  h.  The  Chezy  eoeffieient  ean  be  given  in  terms  of  the  Darey 
frietion  faetor  (fo).  Beeause  the  Fanning  frietion  faetor  is  one  fourth  of  the  Darey  frietion  faetor,  i.e., 
the  Chezy  eoeffieient  ean  also  be  given  in  terms  of  the  Fanning  frietion  faetor  as 
(Henderson  1966) 


C,  = 


(6) 


Equating  5  and  6,  substituting  h  for  R,  and  rearranging  yields  an  equation  (in  metrie  units)  for  fp  in 
terms  of  Manning’s  eoeffieient  and  flow  depth  h  in  meters,  i.e.. 


_2gn 

Jf  ^1/3 


(V) 


The  validity  of  Equation  7  has  not  been  proven,  and  it  is  based  on  the  assumption  that  frietion  faetors 
and  Manning’s  n  assoeiated  with  steady  supereritical  overflow  that  has  reaehed  equilibrium  (e.g., 
Chezy  or  Manning  equation)  will  be  the  same  for  unsteady,  rapidly  varying  flows  due  to  wave 
overtopping.  Therefore,  eaution  must  be  exereised  when  applying  Equation  7. 

Hewlett  et  al.  (1987)  suggested  values  of  Manning’s  n  for  grass-eovered  slopes  exposed  to  steady 
supereritieal  overflow.  They  reeommended  n  =  0.03  for  slopes  of  1:10  (tan  a  =  1/10),  deereasing 
linearly  ion  =  0.02  for  slopes  of  1 :3.  They  reeommended  using  n  =  0.02  for  slopes  steeper  than  1 :3. 
This  linear  relationship  ean  be  expressed  by  a  simple  formula: 

n  =  0.0343-0.043  tan  «  (8) 


for  the  range  1/10  <  tan  a  <  1/3.  Yong  and  Stone  (1967)  reeommended  n  =  0.035  for  steep  grass 
slopes.  Table  2  presents  an  example  of  frietion  faetors  ealeulated  using  Equation  7  for  several  flow 
depths  on  a  1:3  slope  (tan  a  =  1/3). 
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Table  2.  Fanning  friction  factors,  fp,  for  1:3  slope. 

Depth,  h  (ft) 

Depth,  h  (m) 

fp  for  n  =  0.02 

fp  for  n  =  0.035 

0.5 

0.152 

0.015 

0.045 

1.0 

0.305 

0.012 

0.036 

2.0 

0.610 

0.010 

0.028 

As  seen  in  Table  2,  the  value  of  frietion  faetor  is  sensitive  to  the  ehoiee  of  Manning’s  n.  This  is 
unfortunate  beeause  the  frietion  faetor  is  quite  influential  for  determining  flow  veloeity  magnitudes 
aeross  the  levee  erown  and  down  the  proteeted-side  slope. 

FLOW  PARAMETERS  ON  THE  PROTECTED-SIDE  LEVEE  SLOPE:  Both  European 
investigators  derived  theoretieal  expressions  for  the  wave  front  depth-averaged,  slope-parallel  flow 
veloeity  down  the  proteeted-side  slope  based  on  simplifioation  of  the  momentum  equation. 
Sohiittrumpf  and  Oumeraei  (2005)  presented  an  iterative  solution,  whereas  van  Gent  (2002)  derived 
an  explieit  formula.  A  eomparison  between  the  two  solutions  reveals  only  small  differenees  in  the 
result,  and  both  formulations  approaeh  the  same  equation  in  the  limit  as  distanee  down  the  slope 
beeomes  large  (Sehiittrumpf  and  van  Gent  2003).  For  ease  of  applieation,  van  Gent’s  formula  is 
preferred,  and  it  was  given  as 


K 


(9) 


with 


K2  =  (g  sin 


(10) 


^3- 


ff 


2  [h 


82%  '^B2%, 


1/3 


(11) 


(12) 

A3 

and  Sb  is  the  distanee  down  the  slope  from  the  erest  edge,  a  is  the  angle  of  the  proteeted-side  slope, 
and  hBi%  and  ub2%  are  the  flow  depth  and  flow  veloeity,  respeetively,  at  the  proteeted-side  erest  edge. 
For  long  distanees  down  slope,  the  exponential  term  in  Equation  9  vanishes,  and  the  veloeity 
equation  reduees  to 
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^sb2% 

Flow  thickness  perpendicular  to  the  slope  at  any  point  down  the  protected-side  slope  is  found  from 
the  continuity  equation  as 


L 


(13) 


h 

'‘'sb2% 


^82%  ■  ^B2% 


'^sB2% 


(14) 


Equations  1-14  can  be  used  to  estimate  the  wave  overtopping  peak  velocity  and  associated  flow 
depth  over  a  levee  that  is  exceeded  by  only  2  percent  of  the  incoming  waves.  The  main  input 
parameters  are  the  wave  height  (Hmo),  wave  period  (Tp),  and  flood-side  slope  (tan  0)  needed  to 
calculate  the  2-percent  runup  elevation  {RuivX  the  freeboard  (i?c);  the  protected-side  slope  (tan  a); 
the  crest  width  (5);  and  the  Fanning  friction  factor  {fp).  The  following  example  illustrates  application 
of  the  flow  parameter  estimation  technique. 


Example  1:  Irregular  Wave  2%  Overtopping  Flow  Parameters 


Find:  The  overtopping  flow  parameters  exceeded  by  2  percent  of  the  waves  associated  with  an 
average  wave  overtopping  discharge  rate  of  qw  =  0.1  ft^/sec  per  foot  of  levee.  The  levee  is  covered 
with  good-quality  grass  on  the  crest  and  protected-side  slope. 


Given: 

Hfn{) 

6  ft 

=  zeroth-moment  significant  wave  height  [=  Hs\ 

Tp  = 

8  sec 

=  wave  period  associated  with  the  spectral  peak 

tan  0  = 

% 

=  flood-side  levee  slope 

tan  a  = 

1/3 

=  protected-side  levee  slope 

B  = 

10  ft 

=  levee  crest  width 

g  = 

32.2  ft/sec^ 

=  acceleration  due  to  gravity 

fp  = 

0.015 

=  Fanning  friction  factor  (from  Table  2  for  n  =  0.02 
and  h  =  0.5  ft) 

Calculate  Wave  Runup  and  Levee  Freeboard.  Runup  and  levee  freeboard  are  the  key 

parameters  needed  to  start  the  calculation  procedure.  These  parameters  can  be  determined  using 
established  equations  from  Part  VI  of  the  Coastal  Engineering  Manual  (Burcharth  and  Hughes 
2002). 
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First,  calculate  the  Iribarren  number  (surf-similarity  parameter)  for  the  given  wave  eondition  and 
levee  flood-side  slope  (CEM  Equation  VI-5 -2).  The  deepwater  wave  length  based  on  peak  speetral 
wave  period  is 


L  =^^2  (8  see)' =328  ft 

2n  ”  In  ^  ' 


and  the  eorresponding  Iribarren  number  is 


C  = 


ianO 


0.25 


V(6  ft) /(328  ft) 


=  1.85 


(15) 


(16) 


Below  are  the  runup  equations  given  by  CEM  Equation  VI-5-6  (de  Waal  and  van  der  Meer  1992) 
and  the  equations  for  freeboard  obtained  by  inverting  the  average  wave  overtopping  formulas  given 
by  CEM  Equations  VI-5-24  and  VI-5-25  (van  der  Meer  and  Janssen  1995).  Different  equations  are 
used  aeeording  to  the  value  of  Iribarren  number. 

For  <  2: 


^u2%  —1-5 


(17) 


R=- 


^ mO  op 

5.2 


In 


q„  Vtan6'  I 


0.06 


•(/.  n  Yh  Yp) 


(18) 


For  ^op>2: 


^u2% 


(19) 


R  =-^ln 

2.6 


\Yr  Yh  Yh  Yp) 


(20) 


The  “gamma  faetors”  in  Equations  18  and  20  aeeount  for  slope  roughness,  berm  effeet,  shallow 
depth,  and  wave  direetion.  Details  are  in  van  der  Meer  and  Janssen  (1995)  or  the  Coastal 
Engineering  Manual. 

Using  the  set  of  equations  for  h^op  <  2  and  setting  all  the  gammas  equal  to  unity  yields 

Ki;.  =  1-5  4  ■■ff.o  =  1-5  (1.85)(6  fl)  =  16.6  ft 
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R=- 


^ mO  ^op 

5.2 


In 


Vtan^  1 


■(r,  n  7t  7f) 


"  5.2 


(O.lftVsec)  1 


(32.2  ft /sec")  (6  ft) 


3  1.85  0.06 


•(l.0-1.0-1.0-1.0)  =  ll.lft 


Parameters  at  Flood-Side  Levee  Crest  Edge.  Substituting  the  recommended  values  of 
CAh2%  =  0.22  and  Cau2%  =  1  -55  into  Equations  1  and  2,  respectively,  rearranging,  and  noting  that  Hj 
cancels  out  in  both  equations  gives  equations  for  the  flow  thickness  and  flow  velocity  magnitudes 
exceeded  by  2  percent  of  the  incident  waves.  Using  the  above  values  for  Ru2%  and  Rc  yields 

/i,j,=0.22(^„,-A,)  =  0.22(16.6ft-ll.lft)  =  1.21ft 


and 


=  20.6  ft/ sec 


Parameters  at  the  Protected-Side  Levee  Crest  Edge.  The  flow  depth  and  velocity  at  the 
protected-side  crest  edge  exceeded  by  2  percent  of  the  incident  waves  are  found  using  Equations  3 
and  4.  Substituting  the  recommended  value  of  C3  =  0.75  into  Equation  3,  recognizing  that  the 
distance  Xc  =  B  at  the  protected-side  crest  edge,  and  using  the  flow  parameters  already  estimated  at 
the  flood-side  edge  yields  the  estimates 


^-0.75- 
B 


V 


(1.21  ft)exp(-0.75)  =  0.57  ft 


^Bl%  ^A2% 


Bfp 
2h 


'B2%  J 


f 

(20.6  ft /sec) exp 

V 


(I0ft)(0.015)^ 
2(0.57  ft)  ^ 


=  18.1  ft/  sec 


The  friction  of  the  crest  surface  resulted  in  a  velocity  reduction  of  just  over  12  percent. 

Parameters  on  the  Protected-Side  Levee  Slope.  The  flow  velocity  and  thickness  exceeded 
by  2  percent  of  the  incident  waves  varies  with  distance  down  the  protected-side  slope  until  terminal 
velocity  is  reached  (the  supercritical  flow  condition  represented  by  the  Chezy  and  Manning 
equations).  Eirst,  determine  the  values  ofK2,  K3,  and7f4  from  Equations  10, 1 1,  and  12,  respectively. 
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=  (g  sin  =  |(32.2  ft/sec^  jsin  tan '(1/3)  |  =2.17(ft/sec^^ 


^3- 


ff  1 


ni/3 


0.015 


ni/3 


2  (0.57  ft) -(18. 1ft /sec) 


=  0.0899  (sec/ 


fC  2.17(ft/sec^) 

=Ug^y - ^  =  (18.1  ft /sec) - = -6.038  ft /sec 

^3  0.0899  (sec/ ft") 


Now  the  flow  velocity  exceeded  by  2  percent  of  the  waves  can  be  found  at  any  location  down  the 
slope  using  Equation  9.  For  example,  at  20  ft  from  the  protected-side  edge  of  the  crest  (i.e.,  Sb 
=  20  ft)  the  slope-parallel  velocity  is  estimated  to  be 

(-3  ^2  •  ^3  •  h  ) 

2.17  (ft/sec" 

Woo/  1  /-5 

0.0899  (sec/ft") 

-3 12. 17  (ft/sec"  f  I  •  |o.0899  (sec/ft"  •  (20  ft) 

U2%  =  22.0  ft/sec 

"(si,=20) 


+(-6.038  ft/sec)  exp 


And,  the  corresponding  flow  depth  is  given  by  Equation  14  as 


‘■2%, 


\sh=20) 


^B2%  ■  ^B2% 

0.57  ft -18.1  ft/sec 

22.0  ft/sec 

=  0.47  ft 


Finally,  the  terminal  velocity  reached  if  the  slope  is  long  enough  is  estimated  using  Equation  13,  i.e.. 


’■2%, 


K, 


2.17(ft/sec")‘'' 
0.0899  (sec/ft" )‘'' 


24.1  ft/sec 


(Note  in  the  above  calculations  the  importance  of  including  dimensional  units  to  make  sure  the  units 
balance  in  the  final  answer.) 
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The  flow  depths  estimated  at  the  proteeted-side  levee  erest  edge  and  at  a  loeation  20  ft  down  the 
proteeted-side  slope  braeket  the  0.5-ft  depth  assumed  when  the  frietion  faetor  was  seleeted  from 
Table  2.  If  the  depths  are  substantially  different  from  the  depth  assumed  for  the  frietion  faetor,  it 
would  be  prudent  to  repeat  the  eomputation  using  a  revised  friction  factor. 

Table  3  contains  the  results  for  the  above  calculations  along  with  similar  results  for  average  wave 
overtopping  rates  of  =  0.2  and  0.5  ft  /sec  per  foot  of  levee.  The  table  also  gives  estimates  for 
calculations  using  the  value  of  friction  factor  derived  assuming  n  =  0.035.  Doubling  the  average 
wave  overtopping  rate  from  0. 1  ft  /sec  per  foot  to  0.2  ft  /sec  per  foot  increased  the  terminal  velocity 
on  the  proteeted-side  slope  by  only  13  percent.  However,  the  increased  friction  factor  reduced  the 
terminal  velocity  on  the  proteeted-side  slope  by  about  one  third.  Thus,  the  value  of  the  friction  factor 
is  very  influential  in  this  estimation  procedure. 


Table  3.  Additional  parameter  estimates  for  Example  Problem  1. 

<?» 

(ft'/s) 

Rc 

(ft) 

^42% 

(ft) 

^A2% 

(ft/sec) 

f^B2% 

(ft) 

Ws2% 

(ft/sec) 

f^2%(s/j=20) 

(ft) 

^2Vo(sb=20) 

(ft/sec) 

W2%(s6=«) 

(ft/sec) 

fF=0.015 

0.10 

11.1 

1.21 

20.6 

0.57 

18.1 

0.47 

22.0 

24.1 

0.20 

9.7 

1.54 

23.2 

0.73 

21.0 

0.62 

24.6 

27.4 

0.50 

7.7 

1.97 

26.3 

0.93 

24.3 

0.82 

27.5 

31.3 

fp  =  0.045 

0.10 

11.1 

1.21 

20.6 

0.57 

13.9 

0.52 

15.2 

15.3 

0.20 

9.7 

1.54 

23.2 

0.73 

17.0 

0.70 

17.6 

17.8 

0.50 

7.7 

1.97 

26.3 

0.93 

20.6 

0.93 

20.6 

20.6 

Figure  3  plots  the  velocity  of  the  2-percent  runup  wave  leading  edge  as  it  progresses  down  the 
proteeted-side  slope.  Curves  are  shown  for  the  three  values  of  average  wave  overtopping.  The 
horizontal  lines  are  the  estimated  terminal  velocities  (last  column  in  Table  3)  associated  with  each 
overtopping  rate.  All  curves  are  for  the  wave  and  levee  parameters  given  for  the  example  problem. 
Depending  on  the  levee  cross  section,  the  flow  down  the  proteeted-side  slope  may  not  reach  terminal 
velocity  before  the  slope  transitions  into  a  flatter  berm. 
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Figure  3.  Wave  overtopping  flow  velocity  on  protected-side  slope. 


REMARKS:  The  equations  in  this  Technical  Note  primarily  are  used  to  solve  for  the  velocity  and 
flow  depth  peaks  exceeded  by  only  2  percent  of  the  incident  waves.  The  levee  surface  is  subj  ected  to 
the  peak  velocities  only  momentarily,  with  lower  velocities  for  the  rest  of  the  wave  passage.  Thus, 
duration  of  maximum  flow  is  fleeting,  and  erosion  might  be  relatively  minor  for  each  wave  unless 
the  erosion  velocity  threshold  is  quite  a  bit  lower  than  the  peak  velocity.  The  relationship  between 
average  wave  overtopping  rates  and  corresponding  soil  erosion  rates  has  not  yet  been  established. 
Finally,  it  is  apparent  that  a  better  understanding  is  needed  for  specifying  an  appropriate  value  for  the 
friction  factor  for  various  slope  surfaces.  The  methodology  is  sensitive  to  the  Fanning  friction  factor; 
therefore,  until  better  guidance  becomes  available,  it  is  advisable  to  use  lower  estimates  of fp  in  the 
calculations. 
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Example  2:  Irregular  Wave  50%  Overtopping  Flow  Parameters 


Find:  The  overtopping  flow  parameters  exceeded  by  50  percent  of  the  waves  associated  with  an 
average  wave  overtopping  discharge  rate  of  qw  =  0.1  ft^/sec  per  foot  of  levee.  The  levee  is  covered 
with  good-quality  grass  on  the  crest  and  protected-side  slope,  and  the  wave  parameters  are  the  same 
as  Example  1,  i.e., 


Given: 

HmO 

6  ft 

=  Zeroth-moment  significant  wave  height  [=  77^] 

Tp  = 

8  sec 

=  Wave  period  associated  with  the  spectral  peak 

tan  0  = 

1/4 

=  Elood-side  levee  slope 

tan  a  = 

1/3 

=  Protected-side  levee  slope 

B  = 

10  ft 

=  Eevee  crest  width 

g  = 

32.2  ft/sec^ 

=  Gravitational  acceleration 

If  = 

0.015 

=  Fanning  friction  factor  (from  Table  2  for  n  =  0.02 
and  h  =  0.5  ft) 

Calculate  Wave  Runup  and  Levee  Freeboard.  Runup  and  levee  freeboard  will  be  the  same  as 
calculated  in  Example  1,  i.e.,  =  16.6  ft  andi?c  =  H-l  ft. 

Parameters  at  the  Flood-Side  Levee  Crest  Edge.  Substituting  the  recommended  values  for 
mean  flow  parameters  of  CahsoVo  =  0.17  and  CausWo  =  0.94  into  Equations  1  and  2,  respectively, 
rearranging,  and  noting  that  Hs  cancels  out  in  both  equations  gives  equations  for  the  flow  thickness 
and  flow  velocity  magnitudes  exceeded  by  50  percent  of  the  incident  waves.  Using  the  above  values 
for  RuiVo  and  Rc  yields 

*,„,=0.17(A,„,,-^,)  =  0.17(16.6ft-ll.lft)  =  0.94ft 


and 


«,»./.  =  ^/7(l)■94)  =  ^(32.2  ft/sec’)  (0.94)  ^(16.6ft-ll.lft)  =12.5  ft/sec 

Parameters  at  the  Protected-Side  Levee  Crest  Edge.  The  flow  depth  and  velocity  at  the 
protected-side  crest  edge  exceeded  by  50  percent  of  the  incident  waves  are  found  using  Equations  3 
and  4.  Substituting  the  recommended  value  of  C3  =  0.75  into  Equation  3,  recognizing  that  the 
distance  Xc  =  B  at  the  protected-side  crest  edge,  and  using  the  flow  parameters  already  estimated  at 
the  flood-side  edge  yields  the  estimates 
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^550% 


-0.75- 

B 


(0.94  ft)exp(-0.75)  =  0.44  ft 


and 


^B50% 


f 

~  ^A50% 

V 


Bf, 

2h 


'550%  J 


f 

(12.5  ft/sec)  exp 

V 


(I0ft)(0.015)^ 
2(0.44  ft)  ^ 


10.5  ft/sec 


Parameters  on  the  Protected-Side  Levee  Slope:  The  flow  veloeity  and  thiekness  exeeeded 
by  50  pereent  of  the  ineident  waves  varies  with  distanee  down  the  proteeted-side  slope  until  terminal 
veloeity  is  reaehed  (the  supereritieal  flow  eondition  represented  by  the  Chezy  and  Manning 
equations).  First,  determine  the  values  ofK2,  Kj,  and7f4  from  Equations  10, 11,  and  12,  respeetively. 

=(gsina)'^^  =  |(32.2  ft/seo^^sin[tan^' (l/3)]|  =2.17(ft/seo^  j 


K 


3 


'u  1 

1/3 

0.015  1 

7  '  ^B50%  )  _ 

2  (0.44  ft)  .(10.5  ft/sec) 

=  0.1175(sec/ft^)‘^^ 


K  ,  ,  2.17(ft/seeM 

K4  =  Usso%  —^  =  (10.5  ft/see) - - - ^  =  -7.968  ft/see 

^3  0.1175(see/ft") 


Now,  the  flow  veloeity  exeeeded  by  50  pereent  of  the  waves  ean  be  found  at  any  loeation  down  the 
slope  using  Equation  9.  Eor  example,  at  20  ft  from  the  proteeted-side  edge  of  the  erest  (i.e., 
Sb  =  20  ft),  the  slope-parallel  veloeity  is  estimated  to  be 

“50%, ^-^  +  ^4  exp  (-3  ■  Kl  ■  s, ) 

^  -^3 

2.17(ft/sec")''' 

^  0.1175(sec/ft^)" 

-3 (2. 1 7 (ft/sec" )  •  |o.  1 175  (sec/ft" •  (20  ft) 

=17.1  ft/sec 


+  (-7.968  ft/see) exp 


And,  the  eorresponding  flow  depth  is  given  by  Equation  14  as 
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^S50%  ■  ^S50% 

0.44  ft  -10.5  ft/sec 

“50%, 

17.1  ft/sec 

=  0.27  ft 


Finally,  the  terminal  velocity  reached  if  the  slope  is  long  enough  is  estimated  using  Equation  13,  i.e.. 


^50%  “ 

A, 


2. 17  (ft/sec' 

0. 11 75  (sec/ft' y'' 


18.5  ft/sec 


SUMMARY:  This  CHETN  has  summarized  European  research  and  technical  publications 
describing  a  method  for  estimating  wave  overtopping  flow  velocities  and  flow  thicknesses  on  an 
earthen  levee.  The  empirical  equations  are  based  on  small-  and  large-scale  laboratory  experiments, 
and  application  is  limited  to  irregular  wave  overtopping  of  trapezoidal  levee  cross  sections  when  the 
still  water  level  elevation  (i.e.,  surge  elevation)  is  lower  than  the  elevation  of  the  levee  crest.  The 
incident  wave  parameters  are  used  to  estimate  the  elevation  of  the  2-percent  runup  (hypothetical 
runup  elevation  if  the  flood-side  levee  slope  continued  indefinitely).  Once  the  2-percent  runup  value 
is  known,  flow  velocity  and  flow  thickness  can  be  estimated  over  the  levee  crest  and  down  the 
protected-side  slope  as  a  function  of  levee  geometry  (slopes  and  crest  width)  and  an  appropriate 
friction  factor.  Resulting  maximum  velocities  and  flow  thicknesses  are  those  values  that  will  be 
exceeded  by  only  2  percent  of  the  waves.  Different  coefficients  are  used  to  estimate  the  flow  depth 
and  velocity  exceeded  by  50  percent  of  the  waves.  Two  example  problems  illustrate  application. 
These  estimates  are  intended  for  use  in  preliminary  design. 

POINTS  OF  CONTACT :  This  CHETN  is  a  product  of  the  Coastal  Structures  Asset  Management 
Work  Unit  of  the  Coastal  Inlets  Research  Program  (CIRP)  being  conducted  at  the  U.S.  Army 
Engineer  Research  and  Development  Center,  Coastal  and  Hydraulics  Eaboratory  (CHE).  Questions 
about  this  technical  note  can  be  addressed  to  the  author.  Dr.  Steven  A.  Hughes  (Voice:  601-634- 
2026,  Eax:  601-634-3433,  Steven.A.Hughes@erdc.usace. army. mil).  Eor  information  about  the 
research  program,  please  contact  the  Program  Manager,  Dr.  Nicholas  C.  Kraus  (601-634-2016, 
Nicholas. C.Kraus@erdc.usace.army.miP).  Additional  support  was  provided  by  the  Department  of 
Homeland  Security  Eevee  Strengthening  and  Damage  Mitigation  Program.  Beneficial  reviews  were 
provided  by  Bruce  A.  Ebersole  and  William  C.  Seabergh  of  CHE.  This  document  should  be  cited  as 
follows: 

Hughes,  S.  A.  2007.  Estimation  of  overtopping  flow  velocities  on  earthen  levees  due 

to  irregular  waves.  CHETN  Technical  Notes  Collection.  ERDC/CHE  CHETN-III- 

77.  Vicksburg,  MS:  U.S.  Army  Engineer  Research  and  Development  Center. 

An  electronic  copy  of  this  CHETN  is  available  from  http://chl.erdc.usace.army.mil/chetn. 
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